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Abstract. For fermions with degenerate single-particle energy levels, the usual relation between the total 
number of particles and the chemical potential [i is only satisfied for a specific number of particles, i.e. 
those leading to closed shells. The treatment of an arbitrary number of fermions requires a modification of 
the chemical potential, similar to the one proposed by Landsberg for Bose-condensed systems. We study 
the implications of the required modification for fermions in a potential, by calculating the ground state 
energy, the free energy, the density, the partition function and the dynamic two-point correlation function. 
It turns out that the modified relation between the fugacity and the number of particles leads to the correct 
ground state energy and density. But for other quantities like the entropy and the two-point correlation 
functions, an additional correction is required and derived. These calculations indicate that many-body 
perturbation theories based on H — fj,N with Lagrange multiplier \i, are not applicable in unmodified form 
for a fixed number of fermions at low temperature. 

PACS. 05.30.Fk Fermion systems and electron gas 

1 Introduction 

When studying many-body effects in systems with Hamil- 
tonian H for a finite number N of particles, it is often 
convenient to introduce the chemical potential /iasa La- 
grange multiplier associated with TV PT|2||^] , and to inves- 
tigate properties related to H — /iN. E.g., the Gibbs free 
energy G is obtained from e~P G — Tre -,3( ^~ pAr ), where 
the trace is taken over all the degrees of freedom, including 
the number of particles. After the calculations, the chem- 
ical potential is then determined from the constraint that 
the number of particles under consideration be N, e.g., to 
obtain the free energy F from e~" F = Tre~° , where the 
latter trace is taken over the degrees of freedom of the N 
particles, with N fixed. This procedure can provide a very 
efficient way to deal with the finite particle-number, but 
it has to be handled with care |3], because in unmodified 
form it is only valid for closed shell systems. 

We show that a modified expression for the fugacity u 
of the form u = Xe^ 13 , as introduced by Landsberg |Sj in 
the case of Bose condensed systems, instead of the usual 
expression e^, provides a relation between the total num- 
ber of particles and the fugacity for any number of parti- 
cles, also if the ground state is degenerate. It leads to the 
correct ground state energy and density. But even the use 
of the modified constraint on the number of particles does 
not provide correct results for other quantities. E.g. the 
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calculations of the entropy and of the two-point correla- 
tion functions require specific and well-defined corrections, 
to be derived below. 



Our derivation is based on the A^-particle density ma- 
trix, properly anti-symmetrized for fermions. With the aid 
of this quantity we can express the partition function, 
the density and the density-density correlation function 
in terms of single-particle density matrices. The method 
allows for the analytical calculation of the partition func- 
tion and the correlation functions in the low temperature 
limit. We show explicitly how the partition function and 
the density can be calculated in order to elucidate the 
same method used in the more involved calculation of 
the density-density correlation function. In essence, the 
method proceeds as follows. The anti-symmetrization in- 
duces a summation over all possible cycles of the permu- 
tation group. In order to avoid the difficult constraint on 
the cycle lengths, arising from the fixed number of par- 
ticles, a generating function is constructed by forming a 
power series in u with the elements An of the sequence 
{Am\N = 1, . . . , oo} as coefficients, where Am is the re- 
sult of the cyclic summation for N particles. This method 
is standard in quantum statistical mechanics |fil7l8l»nTl| . 
and very useful to construct path integrals for identical 
particles (11112113114115) . Replacing u by | it is indicated 
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as the Z-transform: 

oo 1 
n=0 

Once the generating functions are known we have to 
expand them in a power series in it, and extract the coef- 
ficient of the N th power of u to obtain the expression cor- 
responding with a fermion gas of N particles [T2|. Using 
this technique the partition function, the one-point and 
two-point correlation function can be obtained directly 
i.e. without the use of any quantity related to the chemi- 
cal potential. 

It should be noted that we could easily have adver- 
tised this investigation under the common denominator 
of "non-equivalence between ensembles" [17,18.19. 2 0^211 
I22ll23ll24ll25ll2fill27j . But we refrain from doing so for the 
following reason. Starting from a stationary Hamiltonian 
with a countable spectrum, one can construct a density 
of states for a given total energy and a given number of 
particles. The partition function is the Laplace transform 
of this density of states with respect to the energy, and the 
grand canonical partition function is the generating func- 
tion of this partition function with respect to number of 
particles. In this respect the three functions only describe 
the same information in different parameter spaces 

{E,N} ^ {p,N} ^ {/3,u}. 

Laplace— transform Z— transform 

Only if the transforms are approximately evaluated, the 
warranty that they lead to the same results expires |28| . 
This approach is based on the fundamental transforms 
between Q, Z and S, and it indicates that the calculation 
of properties of a system with a finite number of particles 
using the chemical potential - even with the modification 
- is an approximation as soon as the levels of the Fermi- 
system are degenerate. 

In the next section we introduce the basic methods, 
and illustrate them by analyzing the partition function. 
In the subsequent sections, the same methodology is used 
to study the one-point and the two-point correlation func- 
tions. Finally, we will discuss our results and their impli- 
cations for the iV-body problem. 



at zero temperature) can then be calculated. In the case 
of a one-dimensional (ID) oscillator with frequency u) the 
energy spectrum is given by e v — hui [v + |) , with degen- 
eracy g u = 1. Placing N fermions in the system leads to 
the ground state energy E = Tiuj ^2^ = q (v + |) ■ The oc- 
cupancy density p v of a state expressed in terms of its en- 
ergy becomes p v = (e v < ejv) , where is the lowest 
energy of the unoccupied states. The generalized Heavi- 
sidc function used here is defined as 

q , _ f 1 if x is true 
^ ' ~ 1 otherwise. 

However, the non-degenerate case is the exception rather 
than the rule. We might also consider an illustration of 
the case of open shells. E.g., for a two-dimensional (2D) 
oscillator, the number of states with the same energy t v 
is characterized by the degeneracy g v = u+1, whereas for 
the 3D oscillator g v = (y + 1) [y + 2) /2. The set of states 
with the same energy e„ is called a level or a shell, labeled 
by the (nonnegative integer) index v. 

If L denotes the index of the lowest not completely 
occupied shell (i.e., either empty or partially filled), then 
the number Nl of fermions accommodated in completely 

filled shells is Nl — Y^,u=o 9"- ^ n0 * a ^ shells are com- 
pletely filled, there are N — Nl fermions left over to oc- 
cupy the shell with index L. The ground state energy is 
therefore 

L-l 

E (N) = Y,9»e» + {N-N L )e L . (1) 

The energy difference at low temperature between a gas 
containing N particles and one containing (N—l) particles 
is given by: 

E {N)-E (N-l) = e L , (2) 

provided that the gas with N particles is an open shell 
gas. In textbooks this difference is sometimes used as the 
definition of the chemical potential. In the present paper, 
we follow the convention that \x is the parameter that is 
used to satisfy the constraint arising of the finite number 
of particles as discussed in the introduction. 



2 The partition function of an open shell 
fermion gas 

In this section we will elucidate the mathematical meth- 
ods used throughout the paper, by analyzing the parti- 
tion function of an open shell fermion gas as an example. 
It is well known that the Fermi-Dirac distribution can be 
derived from the probabilistic problem of distributing N 
fermions on a given set of states, assuming that a state can 
contain at most one fermion. The resulting distribution is 
easily derived using generating functions [29j . Assume that 
the energy of a state is given by e v . The energy E of the 
ground state (i.e., the total energy of all occupied states 



2.1 From density of states to generating function and 
vice versa 

The local density of states matrix £1 (E, r', r) for a particle 
in a one-body potential can be defined as 

CO g„ 

Q (E, r', r) = 5 ( E ~ z«) E <P»,k (r') < fc (r) , (3) 

v=0 fc=l 

where ip Ul k (r) are ortho-normal wave functions with cor- 
responding eigenvalues e„. The index v accounts for the 
energy levels and the index k runs from 1 to g Ul where g v 
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is the degeneracy of the energy level e u . The diagonal part The summation runs over all permutations P of the par- 

of this matrix is given by tide indices, and the sign factor (-1) P assures that odd 

x permutations contribute with a negative sign. For con- 

i \ s 1 771 \ i / \i2 , .\ venience in the notations, we limit the treatment to the 

J/ (E, r) = > 6 (E — e„) > \w v & (r) , (4) . , . , r™ i • i i i 

^-^ spm-polarized case. I he non-polarized case can be treated 

along the same lines, see e.g. |50). 

and the density of states fl (E) becomes The partition function for this system of N fermions 

is 

„ CXD 

The single-particle propagator in its spectral represen- The integrals over the configuration space can be done 

tation is usually written as using the decomposition of the permutations in cycles 

oo 9u 7 ,.8,9, 10,31 ,32,33,34 . With Mg denoting the number of 

Ki (r', d\r, 0) = e~ f3e " tp„ k (r') & ( r ) > (6) c y c l es 01 length the partition function becomes 



u=Q k=l 



N 



and it is clearly the Laplace transform of /2(£,r',r) Zp (3\N) = ^ 9 [J2 iMe = N ) x ( 12 ) 

K x {v',3\v,Q) = f Q(E,r',r)e-P E dE, (7) ^ ( _ 1)( <-i)Afc 

J 11 A#„IM/ f ^iWJ) ■ 



£=1 



with /3 adjoined to _E, provided the Laplace transform ex- 
ists, i.e., that the energy levels e v are positive. This means introducing a generating function 
that the origin of the energy spectrum is chosen such that 

it does not exceed the lowest eigen-energy eo. The relation °° 

between the single-particle partition function Z\ (3) and Sp (/3,u) = u N Zp {/3\N) , (13) 

the single-particle propagator n=o 



(/3) = J j K x (r',/3|r,0)5(r'-r) dvdv' , 



the constraint J^eLi ^^-l = N can be removed, and the 
summation over Mi can be done 



implies that the single-particle partition function can also / oo , \t \ 

be written as a Laplace transform of the density of states: ~; p u ) = ex p j _ [ u ' % x (gfi) J _ (14) 



OO r. 

Z x 09) = ]T = °( E ) e'^dE, (8) 

n •/ 



y=0 

with the inverse relation given by 



Using the spectral representation JSJ) of Z\ (3) , the sum 
over the cycles can easily be performed. The result is: 



1 r Z F (3, u) = exp f [ Q(E)bx(l +ue~ l3E ) dE ) (15) 

27TZ ./„ 00 



where c is a Bromwich contour. 

The propagator JCooiN distinguishable particles with 

position vectors r 3 (j = l,...,N), all in the same poten- Again as in the case of a single p articlej the two var i a bles 

tial, is the product of their single-particle propagators ^ u have n0 specific physical meaning with respect to the 

N parameters defining the system, i.e., the energy E and the 

K D (f 3- f' 0) = TT Ki (r' /3|r ■ 0) (10) number of particles N. The partition function Z F (J3\N) 

-j--'- ^ ' for TV particles can then be obtained from its generating 

function Sp (3, u) by the inverse transformation 

The vector f denotes the configuration (ri, . . . , rry) T of 1 d N " (3 ) 

the N particles in the system. When fermions in the same Zp [B\N) — " ~~ ' " ' 



N\ du N 



(16) 



u=0 



spin state are involved one has to anti-symmetrize this 
propagator 

Remark that in the derivation of l|13fl u is nothing but an 
^ N / \ intermediate variable. If one substitutes u by l/z, one can 

K F (f, P\f', 0) = — Y (-!) P II Kl \{ Pr h ' ^ r 'r °) ■ state that S \) is the z transform of the sequence 



P 3 = l 



(11) {l,Zp(d\l),...,Z F (0\N),...}. 
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The inverse of the Z transform is given by a contour- 
integration, leading to the following alternative inversion 
formula: 



z F W\n) = 



1 



c 1 



N+l 



(17) 



where C is a closed counterclockwise contour around the 
origin and chosen in such a way that only the origin is a 
singularity. The inversion can also be done using recursion 
formulas [55ll3l)) . Although they are more appropriate for 
bosons, they are occasionally used for fermions |37j . 



2.2 The inversion method and the chemical potential 

Since SV(/3,7) is an analytic function of 7, a circular 
contour with radius u can be considered, and introduc- 
ing 7 = ue 18 the inversion integral can be rewritten as 



z F (P\n) 



1 

2^ 



e [ln3 P (0,ue*»)-Nlnu} e -iNe dd _ (lg) 



Choosing the radius u such that the argument 

[In ^ ((3,ue ie ) -iVlnu] 

of the exponential function reaches an extremum, one ob- 
tains the condition 



N = u-^-\n~ F (3,u) 
du 



f2(E) 



we 



-PE 



1 + ue-P E 



dE. (19) 



Denoting the solution u obtained from this implicit equa- 
tion by u e , the partition function can be rewritten as 



2jt 



Z F {/3\N)=Z e (P\N) V(0)d0, (20) 
Jo 

3 {f3, u e ) 



Ze (P\N) = 



9(9) 



1 S(/3, 



27T ~(0,Ue) 



Ue& ' e~ iNe 



(21) 
(22) 



The free energy T F {P\N) = -jjlnZp (/3\N) therefore 
consists of two contributions 



to use this method also for this type of 
inversion of a Z transform. In this case, the approximation 
turns out to be rather poor, and in general the error esti- 
mate poses a problem. Below we shall show that, although 
the approximation l|24|) is excellent for fermion systems 
with completely filled shells in the low temperature limit, 
it requires important corrections when the shells are not 
completely filled. It should be emphasized that there is 
no approximation involved in the use of u e (or any other 
integration radius u) as long as only the singularity in the 
origin is contained in the integration domain. The approx- 
imation stems from the replacement of the integrand by 
its extremum, in the hope that the contribution of the 
remaining oscillatory contribution to the free energy is 
negligible. 



2.3 The modified fugacity 

The optimal radius u e for inversion can be calculated from 
(I19JI if the energy spectrum of the system is known. After 
integration over E, one obtains the following transcenden- 
tal equation for u e : 



u e e 



-Pe v 



u e e 



-Pe v - 



(25) 



Assuming the usual relation between the chemical poten- 
tial and the fugacity u e = allows in the low tempera- 
ture limit only solutions that consist of closed shells @]. In 
order to have the possibility to incorporate also open shells 
we assume that u e depends on (3 in the low temperature 
limit as follows: 

u e = Xe 0eL , (26) 

where L is the lowest not fully occupied level. Splitting the 
summation over the levels in the equation for N into three 
parts we obtain 



L-l 



(27) 



9L 



1 + A + ^ 9u i + XefX.tr.-**) •■ 

iy—L-\-l 



1 



F F ((3\N) =Fe(P\N)-jj]nU W{9)d9 



withFo (/3|N) 



1 S(u G ,/3) 



(23) 



(24) 



If one would redefine u e as u e = e^, the relation 119|l 
becomes the common expression for the chemical potential 
11 given the total number of particles N. However, this 
expression is a consequence of the mathematical technique 
used for the inversion of the generating function, and it 
leads to an approximation for the partition function. 

It is well known that the steepest descent method of- 
ten gives good results in the case of inverse Laplace trans- 
forms, but several suggestions are found in the literature 



and taking the limit (3 — > 00 we find the following implicit 
expression for the parameter A: 



N = N L + g L 



A 



1 + A' 



(28) 



where Nl is the total number of states in all fully occupied 
levels 

L-l 



(29) 



This number and therefore the level L itself, has to be 
determined from 



N L < N < N L+1 . 



(30) 
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Once Nl is known, one can calculate the value of A, i.e 

TV — Nl TV -N l 



N L +g L ~N Nl+i-N' 



(31) 



Note that this calculation explicitly assumed that TV dif- 
fers from Nl- In the case that TV — Nl, i.e., if the level 
L is empty, one obtains u e = eP eh in the low temperature 
limit. It is interesting to remark that the divergent terms 
in the steepest descent inversion for bosons could be elim- 
inated by Landsberg using a similar modification of the 
fugacity (SJ- 

Using now the generating function E((3,u e ) in the 
limit (3 — > oo with u e — Ae^ eL , one obtains the follow- 
ing zero-order approximation Z e (/3|TV) for the partition 
function 



(TV — Nl) (Nl+i ■ 



N L + 1 -N L 



N) 



N L + 1 -N 



with the fermion ground state energy Eq given by 



L-l 

r 27T , 



(32) 



(33) 



The integral correction L n & (9) d9 to the partition func- 
tion can be derived from 



lim <P(9) 



1 /1 + Ae s 



2tt 



id \ 9L 



1 + A 



-i9(N-N L ) 



(34) 



and hence, using the binomial expansion of (l + Xe i6 ) 9L , 
one obtains 



lim / V{9)d9 

/3^°° JO 



N. 



L+l 



x _ n N l L )(N- V L )> >W) 



(TV L+1 - N) Nl+1 ~ N 



(N l+ i - Nl) 1 
The partition function thus becomes 



lim Z F (j3\N) 

(3— YOO 



9L 

N-N L 



-pE 



(36) 



where the combinatorial factor in front accounts for the 
number of possibilities to occupy the open shell. In order 
to obtain this factor, the modification of the fugacity is 
not sufficient: the inversion procedure has to be carried 
out completely. In other words: the open shell partition 
function of fermions in the low temperature limit has to 
be calculated by an inversion because the contribution of 
the integral correction turns out to be crucial! 

But the exact ground state energy is obtained 
from the approximate expression Q24[l. This means that 
the integral correction to the free energy l|23(l becomes 
negligible in the low temperature limit; it only contributes 
to the degeneracy factor in the partition function, i.e., to 
the entropy. In the next section we will investigate whether 
this coincidence also remains valid for the correlation func- 
tions. 



3 The density of an open shell Fermion gas 

In this section we derive the one-point correlation function 
that describes the density in the ground-state, using the 
exact inversion. Furthermore, we show that the use of the 
modified fugacity recovers the exact result. 



3.1 The density and its generating function 

The generating function for the Fourier transform n q of 
the density of a TV-fermion system has been derived in 
It is given by 



£n q (0, U) = J2 «q (P, N ) Z F (P\N) NU N , 



(37) 



N=0 



and it can easily be expressed in terms of the density of 
states 



gn q (13, U) 

S F ((3,u) 



f E (u,(3)f2(E,r)e~ tcl - r drdE, (38) 



where Je {u, 0) has the form of an occupation function, 
but u is still the auxiliary variable to construct the power 
series 



,-/3E 



(39) 



3.2 The exact inversion 



Inverting i|37|) for n q (/?, TV) by the contour integral tech- 
nique, and filling out Ga^ ((3, z) in the form Ij38(l 



n q 09, TV) 



1 I 

N Z F (0\N) 



5> (q) 2™ 



(40) 



r I - fa ^ ze ~ 0€ " dz 
x lim (h z. F {Ji,z)- t, rr— , 

u^o / yt ' 1 1 + ze-P^ z N+1 

\z\—u 

the contour integral can be transformed into an angular 
integral. One obtains after some effort 



I ( L ~ X N- 
lim n q n v (q) + 

/3-»oo TV N L +1 



TV - Nr. 



N L 



n L (q)J , 

( 41 ) 

with n v (q) is defined as the Fourier transform of the den- 
sity contribution from the single-particle wave functions 
in the level v 



9, „ 

v (q) = 2J / ( 

fe=i J 



r)| 2 e- 1 



r dv. 



(42) 



This result for the density could have been guessed by 
realizing that for an open shell the states belonging to that 
shell must have the same probability of being occupied. 
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3.3 The optimal fugacity 

The generating function for the density in the spectral 
representation Q can be written as 



-fa* 



Hp {[}, u) ~ ^ 1 + ue-P' 



-n u (q) 



(43) 



Repeating the approximate inversion as for the partition 
function from the previous section, using again 



one finds 



L-l 



lim _ , „ r = 2^ (q) + — ~^~ n L (q) ■ (44) 



/3->oo ^ F (/3,U e ) 



N L+1 - N L 



The modified fugacity thus allows to calculate the ground- 
state density exactly, similarly as we found for the ground- 
state energy. 



3.4 A useful relation 

Because the generating function often can be written in 
closed form, the following exact relation 

limn, (ft = ~ to* ^f^t (45) 

/3— oo W ,3-^oo JZp {p,U e ) 

between the results for the ground-state of the two inver- 
sion methods turns out to be of practical numerical impor- 
tance. The density, i.e. the sum of single-particle contri- 
butions over all the occupied states, can be obtained from 
an appropriate limit of a ratio of two generating functions. 
For fermions in a harmonic potential Qn (j3, u) jSp (/?, u) 
is known in closed form and the limit can be used to obtain 
an explicit result for the density, thus avoiding a tedious 
summation over occupied states. 



4 The polarization function of an open shell 
fermion gas 

In this section we investigate whether the two-point cor- 
relation function that describes the density-density corre- 
lations of the fermion system can also be obtained in the 
low temperature limit using a modified chemical poten- 
tial. Unfortunately, this is not the case. The exact inver- 
sion of the generating function and the use of the chemi- 
cal potential give analogous but slightly different results. 
Both calculations use the same techniques as explained in 
the preceding sections. Therefore only an outline of the 
method and a summary of the results will be given. 

The Fourier transform of the two-point correlation func- 
tion is defined as: 



N N 



\j=l 1=1 I 



(46) 



An outline of the derivation of this function can be found 
in |41) . Defining the following generating function for the 
two point correlation function 



G d , (r|/3, u)=J2 <?q',q (t\0, A) Z F (J3\N) A 



2 u N , 



N=0 



a straightforward application of the techniques illustrated 
in our former publications |13II14| leads after a lengthy but 
straightforward calculation to an expression obtained by 
performing the contour integral of the inversion exactly in 
the limit (3 — > oo : 



lim A 2 C q -, q (r|/3,7V) = A (q,q',r) 

A — Al 



(48) 



-M (q,q',r) 



N L+1 -N L 

A - N L N - N L - 1 
N L+1 - N L N L+1 — Nl — 1 
x M (q, q', t) , 



+ 



with 



L-l L-l 



A (q, q', r) = ^2 n - (0) E n »' (<*') ( 49 ) 
i/=0 v'=0 
L-l oo 

+ E E ^"-^AW (q,q'), 

u=0 v'=L 

L-l 

A 1 (q,q',r)=nL(q)^n I/ (q') (50) 

v=0 
L-l 

+ n L (q') ^ n » (l) 

oo 

+ ^e T ^- e -)ML, I ,(q,qO 



L-l 



^V^-^A^x (q,q') 



l/=0 



M (q, q', r) = n L (q) n L (q') - M L , L (q, q') , (51) 



and 



M v y (q,q')=]T$> (q').(52) 

fe=i k'=i 

A„, k; „<, k , (q) = J < fc (r) p„, ifc , (r) e^dv. (53) 

It should be noted that the expression for A 2 (q, q', r) 
implies that this term is independent of Euclidean time r. 

It turns out that the result l|48[) is not obtained with 
the approximate inversion (r|/3,u e ) / Bp (f3,u e ) of 

the generating function. Using again u e = \e^ tL the dom- 
inant terms for f3 — > oo instead become 



G C , i T \Pi u sd) 

lim — r — = A (q, q , r) 

/3-.OC ^p(P,U s d) 



(54) 
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N 



N L 



N L+1 - N L 
N-N L 



M (q, q',r) 

N-N L 



+ N L+1 ~ N L N L+1 - N L 
x A 2 (q,q',r) . 

Although very similar to the exact result (|48|) the combi- 
natorial factors for the not fully occupied level are slightly 
different. The justification of this approximate inversion is 
rather elaborate. E.g., it is not trivial that the contour in- 
tegration around the origin can be extended to the large 
radius u e without introducing extra poles. This point has 
been carefully checked. 

Although very similar to the exact result l|48|l the com- 
binatorial factors in (|55fl for the not fully occupied level 
are not correct. Therefore the similar correspondence as 
found for the partition function and for the density (|45|) . 
is valid for the two-point correlation function in the case 
of closed shells only. Indeed, Ax (q, q', r) and A 2 (q, q', r) 
do not contribute to the polarization function for the case 
of closed shells, i.e., for N = N^. This leads to: 



lim N 2 C a , 

0->ca 



q%q (t\0,N)= lim 



G c„ 



(r\/3,u e ) 



5 (/3,u e 



(55) 



5 Discussion and Conclusions 

Several authors suggested or stated that calculations per- 
formed with the methods developed in the grand canonical 
ensemble cannot be trusted when one deals with a given fi- 

27 . In this paper we have pinpointed the origin of this 
discrepancy for open shell Fermion systems. Our analysis 
is based on an exact inversion of the existing transform 
between the partition function, the one-point correlation 
function and the two point-correlation function for N par- 
ticles and their generating functions. We observed that 
the generating function of the partition function for the 
proper parameter choice becomes the grand canonical par- 
tition function of the same model. This observation allows 
to compare the results for the exact inversion with those 
relying on the chemical potential. It should be stressed 
that a modification of the relation between the fugacity 
and the chemical potential is needed to obtain the correct 
ground state energy, but this modification is not sufficient 
to obtain the low temperature limit of the partition func- 
tion correctly in the case of open shell fermion gases: the 
entropy requires the exact inversion. In summary the mod- 
ification of the fugacity definition for open shells leads to 
the correct expressions for the ground-state energy and 
the ground-state density. But it is not good enough to ob- 
tain the low temperature limit of the partition function 
and the polarization function, and a correction is needed, 
as derived above. This situation is summarized in table 1. 

Note that the differences discussed here have their ori- 
gin in the degeneracy of the lowest not completely oc- 
cupied level. In the absence of this degeneracy, i.e. for 



Table 1. This table summarizes for which quantities and type 
of systems the use of the unmodified chemical potential /x and 
of its modified form leads to correct results. 



Quantity 


/' 


Modified fj, 


Ground state energy 


closed shell 


open shell 


Density 


closed shell 


open shell 


Partition function 


closed shell 


extra correction 


Correlation function 


closed shell 


extra correction 



closed shell systems, the ground-state expressions of both 
inversion schemes are identical. But if there are N — Nl 
fcrmions in the lowest not completely occupied level, a 
correction term proportional to 

N~N L 
N L+1 - N L 

is introduced in the ground-state energy and the density 
by modifying the fugacity. An extra term proportional to 

(N-N L )(N-N L -1) 



(N L+1 - N L )(N L+1 - N L - 1) 

is found for the polarization function. 

The relation between the density of states, the parti- 
tion function and the generating function have been dis- 
cussed by several authors. In classical statistical mechan- 
ics they are attributed to Fowler and Darwin 42,43 . In 
quantum statistical mechanics they are discussed in the 
book by P.T. Landsberg who makes an epistemologi- 
cal analysis of them (see, e.g., ref. |44j). A mathematical 
justification and review of the present understanding can 
be found in 05] • Of course the present findings also oblige 
to scrutinize the use of the chemical potential in other 
calculation schemes like the density functional approach 
or the Gross-Pitaevsky equation. The questions are there: 
"How can one control these methods, clearly based on the 
grand canonical ensemble, on their correctness for a finite 
number of particles?" , and "Is a modification of the chem- 
ical potential in these cases also necessary?" . Our prelimi- 
nary investigations indicate that in the case that the levels 
are degenerate, one should use the modified chemical po- 
tential in order to get the energy and the density correct, 
and to obtain a fair approximation for the two-point cor- 
relation function. 



This work is performed within the framework of the FWO 
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